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Abstract 

The spectral flow of the overlap operator is computed numerically along a particular 
path in gauge field space. The path connects two gauge equivalent configurations 
which differ by a gauge transformation in the non-trivial class of tt^{SU{2)). The 
computation is done with the SU(2) gauge field in the fundamental, the 3/2, and 
the 5/2 representation. The number of eigenvalue pairs that change places along 
this path is established for these three representations and an even-odd pattern 
predicted by Witten is verified. 
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1 Introduction and summary of results 

In 1982 Witten presented an argument that a chiral SU(2) gauge theory 
with a single left-handed doublet is mathematically inconsistent [1] . His obser- 
vation was that the sign of the fermion determinant for this theory cannot be 
defined to satisfy both gauge invariance and smooth gauge field dependence. 
This obstruction is known as Witten's global SU(2) anomaly. Only for an even 
number of doublets is the theory well-defined since in that case the sign does 
not matter. 

Even though a chiral SU(2) gauge theory in the fundamental representa- 
tion is anomalous, many of the higher representations are not. In terms of the 
diagonal generator T3, normalized as usual such that trT 3 2 = 1/2 in the fun- 
damental representation, only those representations are anomalous for which 

2trT 3 2 = odd integer. (1) 

Alternatively, labeling the SU(2) representations by their highest weight j, the 
anomalous representations are given by 

j = 2/ + ^, 1 = 0,1,2.... (2) 

Besides the fundamental representation with j = 1/2, also the representations 
with j = 5/2, 9/2, 13/2 etc. suffer from Witten's anomaly. The representations 
with integer values of j as well as the remaining half-integer representations 
are anomaly free. 

This pattern is somewhat puzzling. Of course, everyone is familiar with 
the difference between integer and half-integer representations of SU(2). So 
one might anticipate the absence of Witten's anomaly for the integer-valued 
representations. But what causes only half of the half-integer representations 
to be anomalous? The distinction between real and pseudo-real representa- 
tions cannot offer an explanation here, since all half-integer representations 
are pseudo-real. 

Let us briefly recall Witten's argument to see where the criterion (1) comes 
from. First of all, the set of SU(2) gauge transformations is not continuously 
connected, but rather falls into two disjoint homotopy classes. This is mathe- 
matically expressed by 

n 4 (SU(2)) = Z 2 , (3) 

and it means that there exist topologically non-trivial gauge transformations 
which cannot be smoothly deformed to the identity. 

Taking a non-trivial gauge transformation g we can define the linear inter- 
polation 

A M (z,t) = (4) 
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Fig. 1. The simplest example of how the eigenvalue flow as a function of the path 
parameter t could look like. One eigenvalue pair crosses zero at t c and changes places 
at the end of the path. 



between an arbitrary gauge field and its gauge transform A 9 ^ = g(A^ + 
d^)g~ l . The t-dependent gauge field (4) defines a smooth path in configuration 
space. 

Next, consider this path as the background gauge field for the massless 
Dirac operator 

D = Tm (^ + A,) . (5) 

D is anti-hermitian and anti-commutes with 75, so the eigenvalues of D are 
purely imaginary and come in complex conjugate pairs. Furthermore, the spec- 
tra at t — and t — 1 are identical because the gauge fields for these t-values 
are gauge equivalent. 

Witten proved, by invoking an Atiyah-Singer index theorem [2] for a cer- 
tain five- dimensional Dirac operator (the fifth coordinate is essentially the 
path parameter t) , that along the path (4) an odd number of eigenvalue pairs 
{A(t), A*(£)} cross zero and change places: 



- K 

t=0 



J 



k = l,...,n, n = odd. (6) 



t=i 



The simplest scenario with only one eigenvalue pair changing places is sketched 
in figure I. Note that the spectrum of D is discrete because Witten considered 
space-time to be a four dimensional sphere and therefore to be compact. 

Having established this result for the spectral flow, Witten showed that a 
chiral SU(2) gauge theory cannot be properly denned. The integration over the 
fermionic fields in the functional integrals for a theory with a single doublet 
of Weyl fermions yields, up to a sign, the square root of the determinant of 
the Dirac operator D. The sign of the square root is ambiguous and needs to 
be defined by hand. As a result of (6), any smooth and local sign definition 
leads to 



y JdetD(A IM (t = 0)) = (-l)" v /det J D(^(t = l)) 



(7) 



Since n is odd gauge invariance is spoiled and the theory is ill-defined. 

So far the discussion was restricted to the fundamental representation of 
SU(2). It turns out that the integer n in (6) depends on the representation. The 
Atiyah-Singer index theorem predicts n to be odd only for the representations 
with 2trT 3 2 being an odd integer. For all other representations n is even and 
gauge invariance can be preserved. 

Witten's anomaly therefore hinges on the the integer n. It is fairly simple 
to show that the representations with integer-valued j are anomaly free, even 
without the Atiyah-Singer index theorem. Since the integer-valued represen- 
tations are real one can show that each eigenvalue of D is twofold degenerate. 
Hence the exponent of the sign in (7) is 2n (or An, 6n etc.) and the sign auto- 
matically vanishes, independently of n. 

The same argument cannot be applied to the half-integer representations, 
which are pseudo-real. The anomaly pattern for the half-integer representa- 
tions is therefore directly tied to n itself. This leads us to the following ques- 
tion: What is the value of n as a function of j that causes only half of the 
half-integer representations to be anomalous? 

At this point it is important to note that Witten uses the criterion (1) 
only to label the anomalous representations. He does not claim that the value 
of 2 tr T 3 2 coincides with the total number n of eigenvalue pairs that change 
places along the path. In view of the numbers for the first three half-integer 
representations, 

j = \: 2trT 3 2 = l, 

j = §: 2trT 2 = 10, (8) 

j = §: 2trT| = 35, 

it is indeed hard to believe that 10 and 35 eigenvalues change places for the 
3/2 and 5/2 representation, respectively. In fact, the mod- two index theorem 
for skew-adjoint elliptic operators [2] used by Witten does not predict any 
details about the spectral flow of the Dirac operator. All it predicts is whether 
n is even or odd, not n itself. 

It is worth mentioning that there is another way to study Witten's anomaly 
by embedding the gauge group SU(2) in a larger group, SU(3) for example 
[3,4]. By doing so one also finds the representations with j given in (1) to be 
anomalous [5,6]. However, also this method lacks any insight concerning the 
spectral flow of the Dirac operator. The value of n remains unknown. 

In this paper we compute n for the 3/2 and 5/2 representations. Our method 
is conceptually fairly simple. We formulate the path (4) on a discrete space- 
time lattice and compute numerically the spectral flow along this path for a 
suitable lattice Dirac operator. Suitable here means that the Dirac operator 
should satisfy the Ginsparg- Wilson relation [7] in order to preserve exact chiral 
symmetry even for finite lattice spacing. The feasibility of this method has 
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already been demonstrated for the fundamental representation [8] (see also 
[9]). There the spectral flow was found to look roughly as in figure 1, where n 
is 1. 

We find n to be 2 and 3 for the 3/2 and 5/2 representation, respectively. 
Hence the results for the lowest three half-integer representations can be sum- 
marized by the simple relation 

n = j + 2 ■ (9) 

It is tempting to propose that (9) holds for all j. After all, (9) is the simplest 
relation one can imagine that explains the odd-looking anomaly pattern in 
terms of the spectral flow. Needless to say, a rigorous proof of (9) for all j is 
beyond the scope of a numerical calculation and needs to be found elsewhere. 

In addition to our results for n we find some unexpected behavior for the 
spectral flow. In order for an eigenvalue pair to change places, it must cross 
zero at least once along the path. The actual spectral flow for the higher 
representations, however, is more complicated. For the 3/2 representation we 
find 4 zero crossings, one eigenvalue pair crosses zero 3 times, a second one 
once. For the 5/2 representation it is even more complicated: In total 9 zero 
crossings are found. The results suggest that one eigenvalue pair crosses 5 
times, another one 3 times and a third one once, even though our numerical 
approach cannot establish this beyond any doubt. 

These results are obtained for one particular path. This is sufficient for 
determining n, since it is path independent and only depends on whether the 
gauge transformation belongs to the non-trivial class of 7^(577(2)) or not. In 
other words, it is a homotopy invariant. In contrast, there is a priori no reason 
why the number of zero crossings should be path independent. It is fairly 
possible that it depends on the specific starting and end point of the path or 
the details of the path itself. 

To check for this the computation of the spectral flow has been repeated 
for different paths. Randomly generated gauge transformations g in the non- 
trivial homotopy class have been used in the definition (4) of the path. For 
all paths studied we have always found 1 zero crossing for the fundamental 
representation and 4 zero crossings for the 3/2 representation. In particular, 
we have never found the naively expected 2 zero crossings for the 3/2 repre- 
sentation. Unfortunately, the same check could not be carried out for the 5/2 
representation for lack of computer resources. 

Our results for the total number k of zero crossings along the linear inter- 
polation (4) can be summarized by 

k = (j + 2 1 ) 2 = n 2 . (10) 

Whether this relation is valid for all representations j and whether it holds for 
all paths in gauge field space is unknown. Again, answering these questions is 
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beyond the scope of a numerical study and this paper. 

Having already summarized our main results, the rest of the paper is orga- 
nized as follows. In section 2 we give a homotopically non-trivial gauge trans- 
formation and construct the lattice analogue of Witten's path (4). We also 
discuss briefly the relevant spectral properties of the overlap operator [11, 12] 
which was used in the numerical computation. After describing some numer- 
ical details the results for the spectral flow are presented in section 3 and 4. 
Concluding remarks are made in section 5 and additional details can be found 
in the appendix. 



2 Preliminaries 



2.1 A topologically non-trivial gauge transformation 

An explicit example of the non-trivial element of 7T4(SU(2)) has been con- 
structed in [10]. The main point in the construction is the observation that the 
group SU(3) is a non-trivial SU(2) bundle over S 5 . Starting from local trivial- 
izations of that bundle on the upper and lower half-sphere of S 5 , the transition 
function on the equator is a homotopically non-trivial mapping from S A into 
SU(2). In terms of the Pauli matrices o~j,j = 1, 2, 3, the final result in [10] can 
be written as 

*( n ) = TTW* + < ?Tf^f>' (U) 

where <7o denotes the two dimensional identity matrix and n is an element of 
S 4 = {n E 1R 5 , n 2 = 1}. The functions £^ are given by 

fo(n) = 2n 4 , 

fi(n) = 2(n n 2 + n x n z ), 

(12) 

£ 2 (n) = 2(n n 3 - nin 2 ), 

= ng + nf-n|-n§. 

From this definition it follows immediately that 



i=i 




£*?=£"? • (is) 



Hence £ = (£i, £ 2 > £3) is an element of S 2 if n = (n , ni, n 2 , n 3 ) is an element of 
S 3 . Using (13) one can easily check that g(n) is indeed an element of SU(2). 

The coordinates n in (11) are elements of S A , used as stereographic coordi- 
nates of M .To introduce a Euclidean space time lattice we need to write g(n) 
in terms of cartesian coordinates x M . This is done by using the well-known 
relations 
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between cartesian and stereographic coordinates. Here x denotes the Euclidean 
distance of the space-time point x M . Note that in this convention the north 
pole 714 = 1 corresponds to x = 00. 

Using the formulae (14) in (11) one finds the corresponding map from M 
into SU(2) as 



q(x) = — A -Un + iV ai * 3KJJj a - , (15) 

yv ; x 4 + 1 Y x + 1 

At the origin g is identical to —a and it becomes the identity in the limit 
x — > 00. Due to the fourth power in x this limit is practically reached for 
moderate values x. 

The mapping g has various symmetry properties. For instance, it is invariant 
under the transformation x^ —> — x M . In addition, there exists a transformation 
T that transforms g into its inverse. This transformation is explicitly given by 



1 



1 , ,^2^(x) 



Xn — > — £3, X 2 — > Xi, 

T: (16) 
Xi -> -x 2 , x 3 -> x . 

This unitary transformation satisfies T 2 = —1 and we find 

g(Tx)=g-\x) (17) 
for the gauge transformation (15). 



2.2 Definition of the path 

Restricting the coordinates x M in (15) to the points of a lattice defines a 
lattice approximation of the gauge transformation g. As long as the lattice 
spacing a is small enough, this approximation will give rise to a spectral flow 
of a lattice Dirac operator equivalent to (6), even if the topological statement 
(3) does not hold for lattice gauge transformations [8]. When working on a 
finite lattice with periodic boundary conditions, one can simply 'cut' g(x) for 
x^ larger than the lattice size L without spoiling the spectral flow, provided 
L is large enough. 

Let us construct the lattice analogue of Witten's path (4). As the starting 
point we take the constant gauge field 

V{x,fi)=exp{ie{(i)a 3 ) (18) 

with some angles 9{fj) which will be specified later. Using V(x,fi) and its 
gauge transform V 9 (x, fi) = g(x)V(x, fi)g(x + afi)" 1 we define the path 
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U t {x, M ) = Nity 1 { (1 - t)V(x, /i) + tW(x, /i) }, 

(19) 

N(t) = det { (1 - t)V{x, n) + tV 9 {x, fx) } . 

The link Ut(x,fi) is well-defined for all values of t as long as V ^ V 9 . The 
normalization factor N(t) ensures that U t (x,/j,) is an element of SU(2). Notice 
that U t is a smooth function of t. 

The path (19) has a very useful symmetry property Taking into account 
N(t) = N(l — t) one finds for the link variables the relation 

U t (x,/i) = g(x)U 1 - t (x, fi)g(x + aft) -1 , 

(20) 

U t {x, M ) = iV(t)" 1 {(1 - t)V{x, /i) + tV*- 1 (x, M )} , 

i.e. U t and U\-t are gauge equivalent. The gauge field U t differs from U t by 
using the inverse gauge transformation g~ l for the gauge transform of V. 

Taking into account (17) and (20) it is not difficult to prove that the spec- 
trum cr of the well-known Wilson-Dirac operator D shows the symmetry prop- 
erty 

a{D{U t )) = a{D{U X -t)) ■ (21) 

It is therefore sufficient to compute the spectral flow along the path for either 
< t < 1/2 or 1/2 < t < 1. This, of course, is beneficial in a numerical 
computation. 

So far we have restricted the definitions to the fundamental representation 
of SU(2). It goes almost without saying that (15) and (19) immediately de- 
fine corresponding gauge transformations and gauge field paths for the higher 
group representations by applying 

g{x) — > Rj [g(x)\ , 
U t (x,fi) — ► Rj[U t {x,n)], 

where Rj denotes the SU(2) matrix for the representation j. The symmetry 
property (21) for the spectrum is valid for all representations. The explicit 
formulae for R 3 / 2 and R5/2 can be found in appendix A. 

2. 3 Spectral properties of the overlap operator 

The overlap operator [11,12] proposed by Neuberger is an explicit example 
of a Dirac operator that satisfies the Ginsparg- Wilson relation [7]. It therefore 
allows the definition of chiral fermions on the lattice which is necessary in the 
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present context .Q 

In terms of the usual Wilson-Dirac operator D w the overlap operator is 
given by 

aD = l- A{A ] AY l/ \ (23) 
A = 1 - aD w . (24) 

The eigenvalues Aj of aD lie on a unit circle around 1 in the complex plane. 
In addition, the eigenvalues come in complex conjugate pairs A,j,A*, because 
the overlap operator satisfies the hermiticity relation 

Dt = T5 £> 75 . (25) 

If uji denotes an eigenvector with eigenvalue Aj, the vector cD, = -750;^ is an 
eigenvector with eigenvalue A*. Since D is a normal operator, its eigenvectors 
can be chosen to form an orthonormal basis for the fermion fields. 

We are interested in the spectral flow of the overlap operator along the 
path defined in (19). The overlap operator D(U t ) = D t admits to choose 
eigenbases that depend smoothly on the path parameter t. Given such a basis 
Vi(t), Witten's result (6) for the spectral flow means that for an odd number 
of eigenvectors the following holds: 

D v ik (0) = \ k v ik (0) D lVih (l) = \* k v ik (l), k = l,...,n r (2Q) 

Let us assume that the spectrum is non-degenerate at the beginning of the 
path. This can always be achieved by choosing the starting gauge field accord- 
ingly, provided that there is no exact symmetry leading to a degeneracy of the 
spectrum. In that case we can conclude 

v i k (°) = ^ > v ih {l) = nu ik , (27) 
«i fc (°) = Ki 1 ) = K>UJ i k , (28) 

with k, k 1 being at most phase factors. In that sense the eigenvectors have 
interchanged along the path. In short, this result might be expressed in terms 
of the imaginary parts of the eigenvalues as 



Im \ ik 



-Im \ ik 



t=Q 



k = l,..., nj . (29) 



t=i 



The integer rij is odd for the anomalous representations with j = 1/2, 5/2, and 
even for j = 3/2. The subscript j has been added to emphasize the dependence 
on the representation j. 

Notice that rij depends only on the representation j and the topology of the 
continuum gauge transformation g, i.e. on the class of 7r 4 (5'?7(2)) the gauge 



1 For an introduction to exact chiral symmetry and chiral fermions on the lattice 
the reader may consult the overviews [13, 15, 14], for example. 
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transformation belongs to. The integer rij does not depend on the particular 
gauge transformation g. It is also invariant under smooth deformations of the 
path (4) because it cannot change continuously. 

In order for an imaginary part to change sign, one would expect it must 
become zero at least once along the path. Even though eigenvector mixing 
sometimes complicates this picture^], let us define 

kj = number of zero crossings (30) 

along the path in the representation j. From (29) one would conclude that 
kj > rij. As we will see later on, kj is indeed larger than its minimal value rij 
for the 3/2 and 5/2 representation. 

One important issue deserves a comment. Due to the presence of the opera- 
tor (A^A)" 1 / 2 smooth t-dependence of the overlap operator is not obvious. And 
indeed, this is only guaranteed if the gauge fields along the path are smooth. 
More precisely, a differentiable t-dependence of D — and consequently for the 
eigenvectors and eigenvalues of D — holds if A^A has no zero mode along the 
path [16]. This, however, can be achieved by choosing the lattice spacing a 
small enough in the discretization of the continuum gauge transformation g. 

The symmetry property (21) also holds for the spectrum of the overlap 
operator. For the eigenvalues as a function of t it implies 

Im Xi(t) = ±Im Ai(l - t) . (31) 

Here the positive sign holds for eigenvalue pairs that do not change sign, while 
the negative sign applies to the rij pairs that satisfy (29). As a consequence, 
the imaginary part of those rij pairs must cross zero at t = 0.5. Turning 
the argument around, any eigenvalue pair that crosses zero at t — 0.5 is one 
of those rij pairs in (29). Hence, in order to determine rij it is sufficient to 
compute the spectral flow for a small neighborhood around t = 0.5. This will 
be essential for our determination of rij. 

2.4 Computing numerically the spectral flow 

Some general remarks concerning the numerical calculations seem to be 
appropriate before discussing the results. For notational convenience we set 
the lattice spacing a = 1 from now on. 

The computation of the spectral flow was done in two steps. First we 
computed the lowest eigenvalues of D^D using the Conjugate Gradient al- 
gorithm [17, 18]. This gives the absolute value |Aj| 2 of the eigenvalues of D. 

Since they lie on a unit circle around one in the complex plane, their ab- 
solute value and imaginary part are related by simple trigonometry, leading 
to 



2 see also appendix B. 
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Im Ai = ±|AjU/l 



4 



(32) 



Recall that both signs occur because the eigenvalues come in complex conju- 
gate pairs. Our goal here is to detect zero crossings of the eigenvalues, which 
means zero crossings of their imaginary parts. By using (32) only we cannot 
unambiguously decide whether an imaginary part crosses zero because the 
information about the sign is lost. 

A way to follow the sign along the path using the eigenvectors of D^D 
is described in [19]. Without going into the details, we state here that the 
imaginary part of \ can be expressed as a scalar product of properly chosen 
eigenvectors of D^D: 



The brackets stand for the usual scalar product for the fermion fields and 
uj +> i, 0J_ > i, denote the chiral eigenvectors of D^D with eigenvalue |Aj| 2 (Note 
that D^D commutes with 75). In addition, the phases of the eigenvectors can 
be fixed by making use of the reality properties of SU(2). As a result, a zero 
crossing of Im A, is accompanied by a sign change for the matrix element in 
(33), which is easily computed numerically. This method was employed in the 
vicinities where the eigenvalues of D^D are close to zero in order to establish 
whether the imaginary parts change sign or not. 

Any practical implementation of the overlap operator needs an approxima- 
tion for the inverse square root of A^A in (23). We used a truncated expansion 
in Chebyshev polynomials [16]. Depending on the lowest eigenvalue of A^A, 
which was computed numerically first along the path, the highest degree in 
the polynomial approximation was chosen up to 180. 

A truncated Chebyshev expansion results in truncation errors for the eigen- 
values as well as their imaginary parts. A second source of uncertainty is the 
numerical error. Both errors, however, can be controlled analytically and strict 
error bounds can be derived for the eigenvalues as well as for their imaginary 
parts [19]. 

Several times we were rather vague about the size of the lattice spacing a 
and the lattice extension L in order to capture the non-trivial topology of the 
continuum gauge transformation g. Some experimentation revealed that an 
8 4 lattice is sufficient for that purpose. All the results presented in the next 
section were obtained on a lattice of this size. 

Finally, the angles #(/i) in the constant gauge field (18) were set to values 
between 0.37T and 0.67T. Different values were chosen for each direction to 
remove as much degeneracy of the spectrum as possible. The reason for having 
used a constant gauge field as the starting point of our path was to produce 
a gap in the spectrum. If we had used the classical vacuum configuration 
the lowest eigenvalues would have been very close to zero along the entire 
path. In that case checking for zero crossings would have required very small 



Im Xi = (u +! i\D\u- s j) . 



(33) 



11 



truncation and numerical errors and would have been computationally too 
time consuming. 



3 Numerical results 

The spectral flow of D^D for SU(2) in the fundamental representation is 
shown in figure 2. Since the computation for the fundamental representation 
is numerically not very demanding, the eigenvalues were computed for the 
entire interval [0, 1]. This serves as a check of the symmetry property (21) of 
the spectrum, which is clearly visible. 

As can be seen, the absolute value for all but one eigenvalue is never zero. 
Only the lowest eigenvalue becomes zero once at t = 0.5. Hence the imaginary 
part of all but the lowest eigenvalue cannot change sign along the path. The 
lowest eigenvalue pair, however, crosses zero at t — 0.5 (figure 3) and we find 

ni/2 = h (34) 
h/2 = 1. 

The results for j = 3/2 are shown in figs. 4 to 7 and can be summarized 
as follows: Two eigenvalue pairs cross zero at t = 0.5 (figure 5), resulting in 
n 3/2 = 2. This is in agreement with the Atiyah-Singer index theorem which 
requires 113/2 to be even. One of these pairs crosses zero again at t — 0.67 
(figures 6 and 7). Due to the symmetry property of the spectrum, this crossing 
also takes place at t — 0.33. In total we find four zero crossings along the path. 
Hence 

" S/2 = 2 ' (35) 
h/2 = 4. 

Figures 8 to 11 show the results for the 5/2 representation. The spectral 
flow of the first four eigenvalues of D^D (fig. 8) has obvious similarities 
with the spectral flow for the 3/2 representation (fig. 4). This time, how- 
ever, there are five different t- values where eigenvalues are close to zero (t = 
0.24, 0.37, 0.5, 0.63, 0.76). At t — 0.5 three pairs cross zero (fig. 9), leading 
to n 5 / 2 = 3. In addition, there is another zero crossing at i ~ 0.756 and two 
zero crossings around t = 0.63 (figs. 10 and 11). Due to the symmetry of the 
spectrum we finally find 

" 5/2 = 3 ' (36) 
h/2 = 9. 

As predicted by the Atiyah-Singer index theorem, n 5 / 2 is an odd integer. 



12 



CO 



0.2 



O 

0.1 



I I I °°nlnO° I | | | 

0.2 0.4 0.6 0.8 1 

t 



Fig. 2. The lowest four eigenvalues of D^D with SU(2) in the fundamental repre- 
sentation. The numerical and truncation errors are smaller than the size of the data 
points. 
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Fig. 3. The imaginary parts of the lowest eigenvalue pair for SU(2) in the funda- 
mental representation. The eigenvalue pair crosses zero at t = 0.5. 




Fig. 4. The lowest five eigenvalues of D*D for the 3/2 representation. The data were 
taken for 0.5 < t < 1 only and than reflected at t = 0.5. 
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Fig. 5. The imaginary parts of the lowest four eigenvalue pairs for j = 3/2. Two of 
the pairs cross zero at t = 0.5. 
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Fig. 6. The imaginary parts of the lowest three eigenvalue pairs for j = 3/2 around 
t = 0.67. One of the pairs crosses zero at t = 0.67. 
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Fig. 7. A magnification of figure 4. The lowest eigenvalue is well separated from the 
higher ones between the two zero crossings at t = 0.5 and t = 0.67. This implies 
that the pair crossing zero at t = 0.67 crosses zero again at t = 0.5. 
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Fig. 8. The lowest four eigenvalues of D^D for the 5/2 representation. The data 
were taken for 0.5 < t < 1 only and than reflected at t = 0.5. 
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Fig. 9. The imaginary parts of the lowest three eigenvalue pairs for j = 5/2 around 
t = 0.5. All three pairs cross zero. The additional triangular data point at t = 0.5 
corresponds to the fourth eigenvalue. It is unequal to zero and therefore the fourth 
pair does not cross zero. 
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Fig. 10. The imaginary parts of the lowest three eigenvalue pairs for j = 5/2 around 
t = 0.76. There is one zero crossing at t ~ 0.756. 
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Fig. 11. The imaginary parts of the lowest three eigenvalue pairs for j = 5/2 around 
t = 0.63. There is one zero crossing at t « 0.623 and a second one at t « 0.627. 



Note that we are unable to rigorously establish whether one pair crosses 
zero five times, a second one three times and a third one only once. The 
eigenvalues are too close together to disentangle them properly (see fig. 11). 

This emphasizes the importance of the symmetry property (31) for the de- 
termination of rij. Without making use of (31) we would need to keep track of 
each eigenvalue pair along the entire path, carefully counting every zero cross- 
ing. For the 5/2 representation this would have been absolutely impossible, at 
least for our particular path. 

To summarize, the numerical results show the validity of e.q. (9) for the 
first three half-integer representations. A proof for all j is, of course, beyond 
the scope of a numerical study like the one employed here. 

The results for kj are slightly surprising. The number of zero crossings as 
a function of j is given by equation (10), which is larger than required by 
the results for rij. Another peculiarity is that one eigenvalue crossing seems 
to pull down some of the higher eigenvalues. The second and third eigenvalue 
at t ~ 0.67 in figure 6, for example, come very close to zero without actually 
crossing it. 

Before dwelling too much on these results one should make sure that they 
are not purely accidental, related only to the particular gauge transformation 
(15). It seems possible that kj can be made smaller by choosing a different 
gauge transformation g or by smoothly deforming the path. 

All the data shown so far have been obtained on an 8 4 lattice. To check 
whether the peculiarities stem from the finite lattice size L or the lattice 
spacing a, additional computations for the 3/2 representation were performed 
on 10 4 and 12 4 lattices, varying the lattice size and the lattice spacing. The 
results are qualitatively the same as the ones shown here. The main difference 
is that the magnitude of the eigenvalues decreases for larger lattice sizes, as 
expected. 
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gauge orbit of U, 




Fig. 12. Using different home-topically non-trivial gauge transformations gi for the 
path (19) amounts in changing its end point. 

4 Is kj a topological invariant? 

Strictly speaking, the answer to this question is no. The number of zero 
crossings surely is a path dependent quantity. By simply going back and forth 
along the path the number of zero crossings is trivially enlarged. By doing 
so, however, one can only increase kj. More precisely, kj is multiplied by a 
factor which does not depend on the representation. Hence this trivial path 
dependence is easily "divided out" by normalizing to the fundamental repre- 
sentation and the question can be raised for the ratio kj/ ki/ 2 - We are, however, 
more interested in the question whether kj can be decreased, from four to two 
in the case of the 3/2 representation, for example. 

As a first step to answer this question it is certainly useful to check whether 
our result for fc 3 / 2 persists for various non-trivial gauge transformations. Using 
different gauge transformations in the definition of the path means a variation 
of the end point, as sketched in figure 12. If the four zero crossings we found 
above are a mere accident, related only to the particular gauge transformation 
g in (15), we should find examples for the spectral flow with only two zero 
crossings. 

Topologically non-trivial gauge transformations can be generated numer- 
ically the following way. First of all one generates randomly a lattice gauge 
transformation g with SU(2) elements g(x) at each lattice point. This rough 
gauge transformation can be made smooth by slowly minimizing the functional 



Here V M denotes the usual forward difference operator. After sufficiently many 
smoothing steps the resulting gauge transformation represents a lattice dis- 
cretization of a continuum gauge transformation which is either topologically 
trivial or non-trivial. The topological content can be checked by computing 
the spectral flow of D for the fundamental representation. If g is topologically 
non-trivial, there is exactly one eigenvalue pair changing places, according to 
the findings in the previous section. 




(37) 
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Fig. 13. The lowest three eigenvalues of D^D for the 3/2 representation, this time 
with a randomly generated gauge transformation g in the definition (19) of the path. 
Four zeros are clearly visible. 

Proceeding this way 10 suitable gauge transformations were generated and 
used to compute the spectral flow for the 3/2 representation. Suitable here 
means that the lowest eigenvalue of Ac A along the path was required to be 
not smaller than 0.01 for a given g. For eigenvalues that small the highest 
polynomial degree in the Chebyshev approximation for (A^ 1 A)~ 1 / 2 has to be 
of the order of 300 to guarantee truncation errors as small as 10~ 4 . This is 
already computationally quite demanding. Note that the randomly generated 
gauge transformations do not lead to (17) and the symmetry property (31) 
for the spectrum no longer holds. This doubles the amount of computer time 
necessary to compute the spectral flow. 

Figure 13 shows a typical example of the results. Plotted are again the 
eigenvalues of D^D, and the plot should be compared with figure 7. As an- 
ticipated, the spectrum is no longer symmetric under reflections at t = 0.5. 
However, four zeros are visible. A direct computation of the imaginary parts 
shows indeed four zero crossings. A subtlety concerning fig. 13 is discussed in 
appendix B. 

Note that the second and third eigenvalue are well separated from the lowest 
one at t ~ 0.33 and 0.67. The " almost-zeros" in figure 7 are therefore just a 
feature of the particular gauge transformation (15). 

Qualitatively the same results have been found for all 10 randomly gener- 
ated gauge transformations. Even though the spectral flows differ in the details 
for each gauge transformation, all cases show four zero crossings. No example 
with only two zero crossings has been found. 

The spectral flow was also computed using different starting points for the 
path. Different constant gauge fields as well as a constant gauge field multiplied 
with a smooth but random gauge field were used. As before, no example with 
only two zero crossings has been found. 

It would have been desirable to perform the same checks for the 5/2 rep- 
resentation to see whether we always find nine zero crossings. Unfortunately, 
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this test could not be carried out with the computer resources available to us. 

The number of zero crossings is remarkably stable under the path variations 
we have performed. However, one should note that we had to restrict ourselves 
to a rather special subset of gauge fields for the starting point of the path. In 
order to generate a gap for the spectrum, we have always chosen our starting 
field to be constant, except for small random fluctuations. It is possible that 
the number of zero crossings along the path changes for completely different 
starting fields. 

But even with its short comings, the numerical study of this section clearly 
demonstrates that the result k$/2 = 4 is not connected to the analytically 
constructed gauge transformation (15) and its particular symmetry property. 
Whether the number of zero crossings is of topological origin and given by 
a relation like the one suggested in (10) cannot be answered by a numerical 
study. A possible strategy to show this analytically might be the embedding 
method introduced in [3]. 



5 Concluding remarks 

The main motivation for this work was to better understand Witten's cri- 
terion (1) for the anomalous representations. Recall that the integer 2trT 3 2 
takes the values 1, 10 and 35 for the 1/2, 3/2 and 5/2 representation, respec- 
tively. With the results presented here we conclude that 2trT 3 2 is neither the 
number of zero crossings nor the number of eigenvalue pairs that change places 
along the path. A connection between rij and 2trT 3 2 is obviously more subtle. 

Witten's global SU(2) anomaly is widely known as an obstruction which 
forbids a chiral SU(2) gauge theory with an odd number of Weyl fermions in 
the fundamental representation. That higher representations exist which are 
anomaly free has not caught as much attention. However, it remains to be 
proven that chiral gauge theories based on these representations have indeed 
no global anomalies when formulated on the lattice. This is not as obvious 
as in the continuum. Additional difficulties are imposed by a non-zero lat- 
tice spacing, a finite volume with some specific boundary conditions and the 
locality requirement for the overlap operator [16,20]. 

The integer-valued representations can be shown to be anomaly free also 
on the lattice [21]. Since these representations are real, the Dirac operator has 
an extra symmetry leading to a twofold degeneracy of the spectrum. This can 
be used to establish the absence of a global anomaly. 

The pseudo-real half-integer representations are more difficult and the same 
argument cannot be applied. However, the lack of such an argument does not 
a priori rule out an "accidental" degeneracy for some half- integer representa- 
tions. Our numerical results explicitly show that there is no degeneracy for 
the 3/2 representation. A proof that a lattice chiral gauge theory based on 
this representation does not suffer from a global anomaly is still missing. 
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Appendix 

A The 3/2 and 5/2 representation of SU(2) 

The group SU(2) in the fundamental representation is the set of 2 x 2 
matrices U which satisfy 

U- 1 = U\ detU = 1. (A.l) 

These two conditions are satisfied by matrices of the form 

" I , Ki| 2 + |u 12 | 2 = 1. (A.2) 




The first row of the matrix determines the matrix U and the matrix elements 
Mn an d u 12 can be used to parameterize higher representations of SU(2). The 
3/2 representation, for instance, is four-dimensional and the matrix R 3 / 2 [U] 
reads [22] 



u 3 u 






"12 


— v3wiiiii 2 


wnOnl 2 - 2|wi 2 | 2 ) 


-Ml2(|Ml2| 2 - 2|Wn| 2 ) 


\fZu\ x u\ 2 


y/3u u (u* 12 ) 2 


u* 12 (\u 12 \ 2 - 2\u n \ 2 ) 


<l(Ki| 2 - 2\u 12 \ 2 ) 


V3(u* u ) 2 u l2 


-K 2 ) 3 


VSu* n (ul 2 ) 2 


-VS(u* u fu* 12 


Ki) 3 



V 

The six-dimensional 5/2 representation R 5 / 2 [U] is given by 



/ u 5 
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y/hu\\u 12 

-V5u 4 n u* 12 Wii(|«ii| 2 -4|«i 2 | 2 ) 

10<(n* 2 ) 2 V2u 2 n u* 12 (3\u 12 \ 2 - 2\u u \ 2 )) 
10<(^ 2 ) 3 v / 2m 11 K 2 ) 2 (3|m 11 | 2 - 2|« 12 | 2 ) 
V5n n K 2 ) 4 (n* 2 ) 3 (|n 12 | 2 -4|n n | 2 ) 



l 2 f V^iKs) 4 
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V2u* llU l 2 (3\ Ull \ 2 ~ 2|u 12 | 2 ) v^KJ 2 ^ 

-V^Ki) 2 Mia(3|ui 2 | 2 - 2| Mll | 2 ) v^0Ki) 3 «f 2 

K 1 ) 3 (| Mll | 2 -4| Ml2 | 2 ) >/5(«Ii) 4 ui2 
-V^Kx) 4 ^ Kx) 5 y 

Consider elements infinitesimally close to the identity, i.e. 

u n = 1 + i y , Ul2 = f + z f • ( A ' 3 ) 

Keeping terms linear in only, one easily finds the SU(2) generators for the 
3/2 and 5/2 representation, respectively. By exponentiating them again, one 
recovers the representation matrices given above. In fact that was the way 
they were constructed. 



B Zero crossings and eigenvector mixing 

Looking at figure 12 one might be tempted to explain the sign changes of 
the lowest two eigenvalues by three zero crossings of one eigenvalue and one 
zero crossing of the other one. It is easy to smoothly connect the data points 
in figure 12 with a pen according to this assumed eigenvalue flow. By doing so, 
however, the two lowest eigenvalues need to cross twice, once around t = 0.51 
and again around t = 0.61. 

But this eigenvalue crossing does not take place. Figure B.l shows the 
spectral flow around t = 0.61. It is obvious that the eigenvalues do not cross 
as a function of t (the numerical error in this plot is much smaller than the size 
of the symbols). Similarly, the two eigenvalues do not cross around t = 0.51. 
Instead, one eigenvalue crosses zero four times while the second one is unequal 
to zero along the entire path. 
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Fig. B.l. The lowest to eigenvalues for t ~ 0.61. the eigenvalues do not cross. 

This, however, does not contradict our previous result for 77,3/2. Naively one 
might conclude from this that no imaginary part switches sign along the path. 
But one should remember that the sign change for the eigenvalues is defined 
in terms of the eigenvectors by eq. (26), and not by the t-dependence of the 
eigenvalues. And (26) is still true with two imaginary parts changing sign, 
even though the lowest eigenvalue crosses zero four times. 

The reason we have not encountered this difficulty before is again the sym- 
metry property (31). This symmetry forces the eigenvalues to cross zero and 
each other at exactly the same value t = 0.5, thus enabling us to extract rij by 
a local measurement of the spectral flow around t = 0.5. Without this sym- 
metry one has to track the eigenvectors along the entire path and use (26), 
which is much more complicated. 

The absence of a strict eigenvalue crossing is nothing new. It can also be 
found in QCD in the context of string breaking [23,24,25,26], for example. 
There one finds that the energy level E s for the string state, defined by a 
static qq pair connected by a flux tube, is smaller that the energy Emm of the 
two-meson state for qq separations r smaller than the string breaking distance 
r&. For r > r& the situation is reversed, the energy of the two-meson state 
is smaller than the energy of the string state. Naively one might expect a 
distance where E s equals Emm- Mixing of the states for r « r&, however, leads 
to a gap between the two eigenvalues of the Hamiltonian. But even though 
the energy levels do not cross, this does not contradict E s < Emm for r < r& 
and E s > Emm for r > r&. This is analogous to what we discussed above for 
the spectral flow of the Dirac operator. 
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